Manipulation of Vortices by Localized Impurities in Bose-Einstein Condensates 
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We consider the manipulation ol Bose-Einstein condensate vortices by optical potentials generated 
by focused laser beams. It is shown that for appropriate choices of the laser strength and width 
it is possible to successfully transport vortices to various positions inside the trap confining the 
condensate atoms. Furthermore, the full bifurcation structure of possible stationary single-charge 
vortex solutions in a harmonic potential with this type of impurity is elucidated. The case when a 
moving vortex is captured by a stationary laser beam is also studied, as well as the possibility of 
dragging the vortex by means of periodic optical lattices. 



I. INTRODUCTION 

Interactions between localized impurities, or pinning 
centers, and flux lines in type-II superconductors have 
long been of interest in condensed matter physics [l[, 
with much recent work focusing on the pinning effects of 
arrays of impurities Q ■ Similar studies of the interactions 
between a vortex array in a rotating Bose-Einstein con- 
densate (EEC) and a co-rotating optical lattice 0] have 
further contributed to the interest in the physics of ma- 
nipulating one array of topological structures with a sec- 
ond array of pinning sites. Depending on the configura- 
tion, depth, and rotation rate of the optical lattice, struc- 
tural changes to the vortex array may be induced, and 
have now been experimentally observed Q. Furthermore, 
combining an optical lattice with a rotating BEG may en- 
able new investigations of other interesting phenomena, 
such as for example, alterations to the superfluid to Mott- 
insulator transition [5| , production of vortex liquids with 
broken translational symmetry Q, and the existence of 
stable vortex molecules and multi-quantum vortices 
Yet despite these significant advances, the interactions 
between a single vortex and a single pinning site within 
a BEG, and the associated vortex dynamics, arc not fully 
understood and many problems remain unexplored. A 
more complete understanding of such basic interactions 
may be important for the further development of many 
ideas and experiments regarding vortex pinning and ma- 
nipulation, even for the case of vortex arrays. Here we 
undertake a theoretical and numerical study that exam- 
ines the possibility of vortex capture and pinning at a 
localized impurity within the BEG, and the possibility of 
vortex manipulation and dragging by a moving impurity. 

Manipulation of coherent nonlinear matter-wave struc- 
tures 0, Q in trapped BECs has indeed received some 
examination p^ . For example, in the case of nega- 
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five scattering length (attractive) BEGs in a quasi-one- 
dimensional (ID) scenario, numerical analysis shows that 
it is possible to pin bright solitons away from the center of 
harmonic trap. More importantly, pinned bright solitons 
may be adiabatically dragged and repositioned within the 
trap by slowly moving an external impurity generated by 
a focused laser beam [ll|. Alternatively, bright solitons 
might be pinned and dragged by the effective local min- 
ima generated by adiabatically moving optical lattices 
and superlattices [H, [ll] . The case of repulsive interac- 
tions has also drawn considerable attention. In the ID 
setting, the effect of localized impurities on dark solitons 
was described in Ref. [l^, by using direct perturbation 
theory [l^ , and later in Ref. IigI. b y the adiabatic pertur- 
bation theory for dark solitons [131 ■ Also, the effects and 
possible manipulation of dark solitons by optical lattices 
have been studied in Rcfs. p^. [l9l. [20| . 

In the present work, we limit our study of vortex- 
impurity interactions and vortex manipulation to the 
case of a positive scattering length (repulsive) pancake- 
shaped BEG that is harmonically trapped. We envision 
a single localized inipurity created by the addition of a 
focused laser beam [8|, which may in principle be tuned 
either above or below the atomic resonance, thereby cre- 
ating a repulsive or attractive potential with blue or red 
detunings, respectively. We concentrate on the dynamics 
of a blue-detuned beam interacting with a single vortex. 

Our manuscript is organized as follows. In the next 
section we describe the physical setup and its mathemat- 
ical model. In Sec. IIIII we study the static scenario of 
vortex pinning by the localized laser beam by describ- 
ing in detail the full bifurcation structure of stationary 
vortex solutions and their stability as a function of the 
laser properties and the pinning position inside a har- 
monic trap. In Sec. IIVI wc study vortex dragging by an 
adiabatically moving impurity. We briefly describe our 
observations also for the case of single vortex manipula- 
tion using an optical lattice, and touch upon the possibil- 
ity of capturing a precessing vortex by a fixed impurity. 
Finally, in Sec. |V]we summarize our results and discuss 
some possible generalizations and open problems. 
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II. SETUP 

111 the context of BECs at nano-Kelvin temperatures, 
mean-field theory can be used to accurately approxi- 
mate the behavior of matter- waves @. The resulting 
mathematical model is a particular form of the nonlin- 
ear Schrodinger equation (NLS) known as the Gross- 
Pitaevskii equation (GPE) [2ll,[22l. The GPE in its full 
dimensional form is as follows: 

in^t = ~^y^ij + g\ij\^i' + vix,y,z,t)^, (1) 

where ipi^, y, z, t) is the wavcfunction describing the con- 
densate, TO is the mass of the condensed atoms, g = 
ATrh?a/m and a is their s-wave scattering length. The 
time-dependent external potential V{x,y,z,t) acting on 
the condensate is taken to be a combination of a static 
harmonic trap (HT) holding the condensed atoms, and a 
localized impurity (Imp) provided by a narrowly focused 
laser beam: 

V{x,y,z,t) = VHT{x,y,z) + Vimp{x,y, z,t). (2) 
Herein we consider a harmonic trap potential 

Vht{x, y, Z) = (^^.2 ^ ^ (3) 

with trapping frequencies lo^ and iOz in the radial and z 
directions respectively. In general, Vimp can be a negative 
or positive quantity, corresponding to an impurity that is 
an attractive or repulsive potential for the trapped atoms. 

In the present study we further limit our atten- 
tion to quasi-two-dimensional condensates, the so-called 
pancake-shaped condensates, by considering that lo^ S> 
LOr and that the tight (z) direction condensate profile is 
described by the harmonic trap ground state in that di- 
rection [1, Q . We also consider only cases where V[mp 
is only a function of x and y, and possibly t, and here- 
after remove the z dependence from our notation. Un- 
der this assumption it is possible to reduce the three- 
dimensional GPE ([T]) to an effective two-dimensional 
equation that has the same form as its three-dimensional 
counterpart but with g replaced by g2D = g/^f^az^ 
where = y/h/TrrujJz) is the transverse harmonic os- 
cillator length 

Furthermore, by measuring, respectively, two- 
dimensional density, length, time, and energy in units of 
^z/ff2D, 0.Z, ^z^, and hjjz, one obtains the standard 
form for the adimensionalized GPE in two dimensions: 

iut = -]^{u.:^x + Uyy) + \u\^u + V{x, y, t)u, (4) 

where the harmonic potential now reads 

VM^,y)^^{x' + y^), (5) 

and CI = ujr/iOz is the adimensionalized harmonic trap 
strength. We use throughout this work a typical value 



for the harmonic trap strength of = 0.065 unless stated 
otherwise. Other harmonic trap strengths gave qualita- 
tively similar results. In addition to the harmonic trap 
we impose a localized potential stemming from an ex- 
ternal localized laser beam centered about {xa(t),ya(t)) 
that in adimcnsional form reads 

Virnp{x,y,t) = V^J^ exp I — 1 . 

(6) 

In this equation, V-^^^^ is proportional to the peak laser 
intensity divided by the detuning of the laser from the 
atomic resonance, and e = wo/V2 where 2wo is the adi- 
mcnsional Gaussian beam width. A positive (negative) 
^mp corresponds to the intensity of a blue-(red-)detuned, 
repulsive (attractive) potential. 

Steady-state solutions of the GPE are obtained by sep- 
arating spatial and temporal dependencies as u{x, y, t) — 
'^{x,y) e~*^*, where ^I* is the steady-state, time- 
independent amplitude of the wavcfunction and ^ is its 
chemical potential (taken here as ^ = 1 in adimcnsional 
units). Under the conditions that the location of the 
impurity is time-independent such that Xa and j/a are 
constant, this leads to the steady-state equation 

= -\ + *,,) + 1*1'* + V{x, y)^. (7) 

The initial condition used in this study was one that 
closely approximates a vortex solution of unit charge s = 
±1 centered at (xo,j/o): 

■^{x,y) ^ -^Tvix^y) tanh[(a; - x^f + {y - ya)^] 

n (8) 

X exp [is tan {{y - yo)/{x ~ xo)}\ , 

where '^Tp{x,y) = y^ma.x{fj, — V{x, y), 0) represents the 
shape of the Thomas-Fermi (TF) cloud formed in the 
presence of the relevant external potentials Subse- 
quently, this approximate initial condition was allowed to 
converge to the numerically "exact" solutions by means 
of fixed point iterations. 

III. THE STATIC PICTURE: VORTEX 
PINNING AND THE BIFURCATIONS BENEATH 

A. Vortex Pinning by the Impurity 

It is well-known that a vortex interacting with a har- 
monic trap undergoes a precession based upon the heal- 
ing length of the vortex and the parameters which de- 
fine the trap [H, M, M, M, M, m,M,^- Since we 
arc introducing a localized impurity into the trap it is 
worthwhile to first observe the behavior of a vortex in- 
teracting with only the localized impurity, in the absence 
of the harmonic potential. We note that in this case, the 
parameters Xa and ya may be neglected and the param- 
eters xq and yo can be interpreted as the coordinates of 
the vortex relative to the impurity. 
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FIG. 1: (Color online) Density plot showing a snapshot of 
the interaction of the vortex with a localized impurity in the 
absence of the harmonic trap (i.e., f2 = 0). The presence of 
the impurity (wider field depression) induces a clockwise ro- 
tation of the vortex (narrower field depression) along a path 
depicted by the dark dots. The parameters are as follows: 



(/i, n, s, 



(0) 



(1, 0, 1, 5, 1). The colorbar shows the con- 



densate density in adimensional units. 




FIG. 2: (Color online) Phase diagram depicting the vortex 
pinning by the localized impurity of strength VjJ^p and width 
e. Each curve represents a different pinning location at the 
indicated radii (i.e., distance from the center of the harmonic 
trap). Parameters are as follows: (/i, fi, s) = (1,0.065,1). 



By symmetry, a vortex placed at the center of an im- 
pm'ity (i.e., xq = yo = 0) will result in a steady state 
without precessing. However, a vortex placed off cen- 
ter with respect to the impurity will precess at constant 
speed around the impurity due to the gradient in the 
background field induced by the impurity [3l| . In order 
to study this behavior in a simple physically meaningful 
setting we start with a positive-charge (s = 1) vortex 



without the impurity and then the impurity is adiabati- 
cally switched on at a prescribed distance away from the 
center of the vortex. We find that for V^^^ > the vortex 
then begins to precess around the impurity in a clockwise 
direction. Reversing the sign of V^^^ in order to create 
an attractive impurity induces a counter-clockwise pre- 
cession with respect to the impurity. An example of the 
vortex precession induced by the impurity is shown in 
Fig. [1] It is crucial to note that if the impurity is turned 
on "close enough" to the steady-state vortex such that 
the impurity is within the vortex funnel then the vortex 
would begin its usual rotation but would be drawn into 
the center of the impurity, effectively pinning the vor- 
tex. This effective attraction is related to the emission 
of sound by the vortex when it is inside the funnel of the 
impurity as described in Ref. [33 |. 

Throughout this work we follow the center of the vor- 
tices by detecting the extrema of the superfluid vorticity 
uj defined as a; = V x Vg where the superfluid velocity in 
dimensional units is [1, S3l 



2 m 



(9) 



where (•)* stands for complex conjugation. 

We now consider the net effect of the pinning induced 
by the impurity and the precession induced by the har- 
monic trap. Since one of our main goals is to find con- 
ditions needed for the manipulation of vortices using 
the repulsive impurity, a minimum requirement would 
be that the impurity's pinning strength is sufficient to 
overcome the precession inside the trap and thus pin the 
vortex very close to the location of the impurity (i.e., 
{xq, uq) ~ {xa, Ua))- Therefore, we seek to find the mini- 
mum conditions that an off-center vortex, at a particular 
radius measured from the center of the harmonic trap, 
could be pinned by a localized impurity at that same 
location. For certain combinations of beam parameters 
(strong beam intensity, or large beam widths), the vor- 
tex will remain localized near this point. For other pa- 
rameters (weak intensity, small beam widths), the beam 
cannot overcome the vortex precession induced by the 
harmonic trap and the vortex would not remain localized 
near the beam position. This would give us a lower bound 
for the possible beam intensities and widths for which a 
vortex might be dragged to the corresponding position 
within the BEC. The existence of such pinned states was 
identified by searching in the impurity parameter space 
(strength Vj|^p and width e) for several off-center radii, 
i.e., distances measured from the center of the trap. The 
results are shown in the phase diagrams of Fig. [2] where 
each curve corresponding to a different radius (decreas- 
ing from top to bottom) depicts the boundary in (Vj^p, e) 
parameter space for which pinning is possible. In other 
words, for the points in parameter space below a given 
curve, one gets primarily vortex precession dynamics in- 
duced by the harmonic trap, whereas above these curves 
(i.e., for strong or wide enough impurities), the vortex is 
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FIG. 3: (Color online) The top panel is a one-dimensional slice 
(for ya ~ 0) of the solution surfaces (represented in the inset 
panel) as a function of the position of the impurity {xa,ya) 
for fixed VJ^p. The vertical axis corresponds to the L^-norm 
squared of the solution (i.e., the normalized number of atoms). 
The dashed lines (yellow (lighter) surface) correspond to un- 
stable solutions, while the solid lines (red (darker) surfaces) 
correspond to stable solutions. The critical value of the ra- 
dius, 3.72 as well as the characteristic value Xa = 2 
and Xa = are represented by vertical lines. The bottom left 
panel shows the squared bifurcating eigenvalues, A^, along 
these branches (A'^ > corresponds to an instability), while 
the bottom right shows the saddle-node bifurcation as a func- 



tion of Vj^-p for Xa = 2 fixed (there are lines at Vj^^ = 0.57 



(cr) 



^imp ~ 0.8, and V^^'^ = 1). The solutions and spectra for each 
of the three branches, represented by circles (and letters) for 
each of the characteristic values of Xa and l^J'^p are presented 
in Fig. |4] For these branches {fi, s, e) = (1, 1, 0.5). 
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trapped by the impurity and stays very close to it. 



B. Steady-state bifurcation structure 

In this section wc elaborate our investigation of the 
pinning statics and the associated dynamical stability 
picture. In particular, wc thoroughly analyze the bifur- 
cation structure of the steady states with single-charge 
vorticity in the setting investigated above (i.e. solutions 
to Eq. ([7])) including their stability. The latter will be 
examined by the eigenvalues of the linearization around 
the steady state. Upon obtaining a steady state solution 
^' of Eq. ([7]) and considering a separable complex valued 
perturbation u — a{x^y)e^* -\- b*{x^y)e^ * of the steady 
state, we arrive at the following eigenvalue problem for 
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FIG. 4: (Color online) The solutions (left, zoomed in to the 
region where the vortices live) and corresponding linearization 
spectra in a neighborhood of the origin (right), for the various 
parameter values indicated by the circles in Fig. |3l The top 
(bottom) row corresponds to V^^^ — 1 (Vj^p = 0.8). 

As evidenced by the bifurcation diagrams presented in 
Fig. [21 there exist three solutions (i.e., steady-state vor- 
tex positions) for any impurity displacement radius in 
the interval (0, a;a^cr), letting j/q = without any loss of 
generality. At the low end of the interval (i.e., for an im- 
purity at the center of the trap, (x^, ya) = (0, 0)) there is 
a transcritical bifurcation where the left (i.e., ones for 
negative Xa) and right (i.e., for positive Xa) solutions 
collide (see below for further explanation) and exchange 
stability. At the other end of the interval (i.e., for an im- 
purity at {xa,ya) = (a^a,cr,0)) there exists a saddle-node 
bifurcation where two new steady-state vortex solutions 
can be thought of as emerging as Xa decreases to values 
Xa < Xa.cr, OT convcrscly can be thought of as disappear- 
ing as Xa increases to values Xa > Xa,cr- Among them, 
one stable vortex position is found very close to the im- 
purity (see cases "A" and "D" in Figs.[n]and Fig.d]), and 
another unstable vortex position further away from the 
trap center (see cases "B" and "E") in Figs.[3]and Fig. [4]). 
Considering only Xa, with ya — 0, two solution branches 
are stable: one with the vortex sitting very close to the 
impurity (see cases "A" and "D" in Figs. [3] and Fig. [5]), 
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and one with the vortex close to the center of the har- 
monic trap (see cases "C" and "F" in Figs. [3] and Fig.d]). 
The other solution is unstable (see cases "B" and "E" in 
Figs. [3] and Fig. 2]), with the vortex sitting in the small 
effective potential minimum on the side of the impurity 
opposite the center of the harmonic trap. This branch of 
solutions collides with the one that has a vortex pinned 
at the impurity and the two disappear in a saddle-node 
bifurcation as the attraction of the impurity becomes too 
weak to "hold" the vortex (this bifurcation corresponds 
to the curves of critical parameter values in Fig. At 
Xa = 0, the potential becomes radially symmetric and 
the solution with the vortex on the outside of the impu- 
rity becomes identical (up to rotation) to the one close 
to the origin. Indeed, at this point there is a single one- 
parameter family of invariant solutions with the vortex 
equidistant from the origin (i.e., for any angle in polar 
coordinates), in addition to the single solution with the 
vortex centered at (0,0). The solutions in this invari- 
ant family, not being radially symmetric and being in a 
radially symmetric trap, necessarily have an additional 
pair of zero eigenvalues in the linearization spectrum to 
account for the additional invariance, i.e. they have 4 in- 
stead of 2 zero eigenvalues (note that the solution with 
the vortex in the center for Xa = only has a single pair of 
zero eigenvalues due to its radial symmetry) . For Xa < 
the branches exchange roles (transcritical bifurcation) as 
the previously imaginary pair of eigenvalues for the sta- 
ble Xa > branch (sec cases "C" and "F" in Figs. [3] and 
Fig. m emerges on the real axis and the pair of previ- 
ously real eigenvalues from the unstable Xa > branch 
(see cases "B" and "E" in Figs. [3] and Fig. |4|) emerges on 
the imaginary axis, i.e. the branches exchange their sta- 
bility properties, becoming the reflected versions of one 
another; see the bottom left panel of Fig. [3] and the right 
column of Fig. [4] In summary, for an impurity that is 
strong enough and close enough to the trap center (specif- 
ically, < Xa < Xa,cr, whcrc Xa,cr dcpcuds of thc Strength 
of the impurity) it is possible to stably pin the vortex very 
close to the impurity. However, if the impurity is too far 
away from the trap's center (specifically, Xa > Xa,cr) the 
vortex can no longer be pinned by the impurity. The top 
image of Fig. [3] depicts the bifurcations via the L^-norm 
squared (i.e., thc normalized number of atoms) of the so- 
lution as a function of the radius, Xa (for ya = 0), and as 
a function of arbitrary impurity location (x^, ya) (inset). 
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FIG. 5: (Color online) Parameter regions for successful ma- 
nipulation of a single vortex inside a harmonic trap. The 
area above each curve corresponds to the successful drag- 
ging region. The bottom panel depicts a zoomed region of 
the top panel where the asterisk and cross correspond, re- 
spectively, to the manipulation success and failure depicted 
in Fig. [6l These panels indicate that higher intensity beams, 
and broader beams, can successfully drag vortices over shorter 
timescales than weaker, narrower beams. 



In fact, the above picture holds for any fixed, suf- 
ficiently large V^^^. Conversely, the same bifurcation 
structure can be represented by a continuation in the 
amplitude of the impurity, V^^^^^ (see bottom right panel 
of Fig. [3]). In particular, for a fixed radius of Xa — 2, 
a continuation was performed and a saddle-node bifur- 
cation appears for V^^^^^ ~ 0.57 where an unstable and 
a stable branch emerge (corresponding exactly to those 
presented in the continuation in Xa)- One can then infer 
that the critical V^^^^ decreases cLS X(i decreases. 



IV. THE DYNAMICAL PICTURE: 
AND CAPTURING 



A. Vortex Dragging 



DRAGGING 



We would like now to take a pinned vortex and adi- 
abatically drag it with the impurity in a manner akin 
to what is has been proposed for bright [Hi [13, [ill and 
dark solitons [3 [3 HOI in the quasi-lD configuration. 
Manipulation of the vortex begins with the focused laser 
beam at thc center of the vortex. The laser is then adi- 



FIG. 6: Successful (top row) and failed (bottom row) vortex dragging cases by a moving laser impurity corresponding to the 
parameters depicted, respectively, by an asterisk and a cross in the bottom panel of Fig. \E\ In both cases, the laser impurity, 
marked by a cross, with (Vj^^p, e) = (5, 1) is moved adiabatically from (0, 0) to (5.43, —5.43) and the snapshots of the density 
are shown every 0.5t*. The top row corresponds to a successful manipulation for r = 18.5 while the bottom row depicts a failed 
dragging for a slightly lower adiabaticity of r = 18. 
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FIG. 7: (Color online) Vortex dragging by an optical lattice potential as in Eq. (|12p . The central high-intensity maximum of 
the optical lattice (depicted by a cross in the panels) is moved adiabatically from the initial position (a;a(0), ya(0)) = (0,0) to 
the final position (5.43,5.43). The top row depicts the BEC density (the colorbar shows the density in adimensional units) 
while the bottom row depicts the phase of the condensate where the vortex position ("plus" symbol in the right column) can 
be clearly inferred from the 2iv phase jump around its core. Observe how the vortex loses its guiding well and jumps to a 
neighboring well in the right column. The left, middle and right columns correspond, respectively, to times t = 0, t — 0.5t* 
and t = 2.5i*. The remaining parameters are as follows: {Vo,k,T,9x,0y) — (1.4,0.3215,20,0,0). 



abatically moved to a desired location while continually 
tracking the position of the vortex, Adiabaticity for the 
motion of the impurity is controlled by the adiabatic- 
ity parameter t controlling the acceleration of the center 



{xa(t),ya{t)) of the impurity as: 



Xa{t) = - -{xi ~ Xf) [I + tanh 



ya{t) = Vt- i^iVi- Vf) [1 + tanh 



t - t* 

T 

t - t* 



, (10) 
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where the initial and final positions of the impurity are, 
respectively, {xi, Ui) and {xf,yf). We will assume ya{i) = 
(i.e., Vi ^ Uf ~ 0) for the discussion below. The instant 
of maximum acceleration is 

t* = tanh"^ (Vl - St^ t, (11) 

where J is a small parameter, S = 0.001, such that the 
initial velocity of the impurity is negligible and that 
Xa(0) ~ Xi and Xa{2t*) w Xf (and the same for y). 
This condition on t* allows for the reduction of param- 
eters and allows us to ensure that we begin with a lo- 
calized impurity very close to the center of the trap [i.e., 
{xa{0),ya{0)) w (0,0)] and that we will drag it adiabat- 
ically to {xf,yf) during the time interval [0,2t*]. The 
next objective is to determine the relation between adi- 
abaticity and the various parameters such as strength 
(^mp) O'' the width (e) of the impurity in order to suc- 
cessfully drag a vortex outward to a specific distance from 
the center of the harmonic trap. In our study we set this 
distance to be half of the radius of the cloud (half of the 
Thomas- Fermi radius). We use the value t* to also de- 
fine when to stop dynamically evolving our system. In 
particular, we opt to continue monitoring the system's 
evolution until tj = 3t*. This choice ensures that a vor- 
tex that might have been lingering close to the impurity 
at earlier times would have either been "swallowed up" 
by the impurity and remain pinned for later times, or will 
have drifted further away due to the precession induced 
by the trap. 

Applying this technique, along with a bisection method 
(successively dividing the parameter step in half and 
changing the sign of the parameter stepping once the 
threshold pinning value is reached) within the span of 
relevant parameters yields the phase diagram depicted in 
Fig. O The various curves in the figure represent the pa- 
rameter boundaries for successful dragging of the vortices 
for different impurity widths (increasing widths from top 
to bottom). All the curves for different widths are quali- 
tatively similar corresponding to higher values of the adi- 
abaticity parameter as the width is decreased. This trend 
continues as s approaches the existence threshold estab- 
lished in Fig. [21 In Fig. [S] we depict snapshots for the two 
cases depicted by an asterisk (successful dragging) and a 
cross (failed dragging) in the lower panel of Fig. [S) 

All of the numerical simulations discussed above deal 
with dragging the vortex by means of the localized im- 
purity. As with previous works of vortex manipula- 
tions we also attempted to produce similar results via 
an optical lattice (OL) potential generated by counter- 
propagating laser beams Q. In one dimension, the case 
of bright solitons manipulated by OLs has been stud- 
ied in Refs. [l^, [l3| while the dark soliton case has been 
treated in Refs. [ll, [H, H^l- For a ID OL, simply de- 
scribed by Vq^{x) = Vbcos^(fca; -I- 9x) where and 9^ 
are the wavenumbcr and phase of the OL, the potential 
minima (or maxima) are isolated from each other pro- 
viding good effective potential minima for pinning and 



dragging. On the other hand, when expanded to 2D, the 
OL reads: 

V^Eix,y) = Vo [cos\kxait) + 9,) + cos\kya{t) + By)] , 

(12) 

where k and 9x^y are, respectively, the wavenumber and 
phase of the OL in the x and y direction. Here we observe 
that each 2D minimum (or maximum) is no longer iso- 
lated, and that between two minima (or maxima) there 
are areas for which the vortex can escape (near the saddle 
points of the potential). This is exactly what we observed 
when attempting to drag a vortex using the 2D OL (fT^ 
without sufficient adiabaticity. The vortex would mean- 
der around the various facets of the lattice outside of our 
control. To overcome this one needs to displace the po- 
tential with a high degree of adiabaticity. In doing so, we 
were successful in dragging the vortices under some re- 
straints (relatively small displacements from the trap cen- 
ter) . An example of a partially successful vortex dragging 
by an OL with potential (fT^ is presented in Fig. [71 As it 
can be observed from the figure, the vortex (whose cen- 
ter is depicted by a "plus" ) is dragged by the OL (whose 
center is depicted by a cross) for some time. However, 
before the OL reaches its final destination, the pinning 
is lost and the vortex jumps to the neighboring OL well 
to the right. This clearly shows that vortex dragging 
with an OL is a delicate issue due to the saddle points of 
the OL that allow the vortex to escape. Nonetheless, for 
sufficient adiabaticity, with a strong enough OL and for 
small displacements from the trap center, it is possible 
to successfully drag the vortex. A more detailed study of 
the parameters that allow for a successful dragging with 
the OL (i.e., relative strength and frequency of the lattice 
and adiabaticity) falls outside of the scope of the present 
manuscript and will be addressed in a future work. 



B. Vortex Capturing 

A natural extension of the above results is to investi- 
gate whether it is possible to capture a vortex that is al- 
ready precessing by an appropriately located and crafted 
impurity. This idea of capturing, paired with the drag- 
ging ability, suggests that a vortex created off-center, 
which is typically the case in an experimental setting, 
can be captured, pinned and dragged to a desired loca- 
tion cither at the center of the trap or at some other 
distance off-center. We now give a few examples demon- 
strating that it is indeed possible for a localized impurity 
to capture a moving vortex. The simulation begins with 
a steady-state solution of a vortex pinned by an impu- 
rity at a prescribed radius and a second impurity on the 
opposite side of the trap at a different radius. Initial 
numerical experiments have been done to determine the 
importance of the difference in these distances from the 
trap center. 

As is shown in Fig. [5] the capturing impurity must be 
located sufficiently lower (i.e., closer to the trap center) 



8 



10.00 -r 



5.00 - 



y 0.00 - -I- 



-5.00 - 




-10.00 I 1 1 \ \ 1 \ 1 

-10.00 -5.00 0.00 5.00 10.0 

X 

FIG. 8: (Color online) Capturing a precessing vortex by a sta- 
tionary impurity. The different paths correspond to isolated 
vortices that are released by adiabatically turning off a pin- 
ning impurity at the following off-center locations: (5, 0) (thin 
black line), (6.5,0) (thick red line) and (8,0) (blue dashed 
line). The capturing impurity is located at (—4,0). The first 
and third cases fail to produce capturing while the second 
case manages to capture the vortex. One interesting feature 
that we observed in the case of successful capture is that, 
before its gets captured, the vortex gets drawn into the im- 
purity potential, but then almost gets knocked back out by 
the phonon radiation waves created from the capture which 
bounce around within the condensate. 



than the trapping impurity in order for the vortex to be 
puUed from its precession and be captured by the impu- 
rity. Intuitively one might come to the conclusion that if 
the vortex and impurity were located the same distance 
away from the center of the trap, then the vortex should 
be captured. But due to the interaction between the vor- 
tex and the impurity that was discussed earlier, as the 
vortex approaches the impurity, it begins to interact with 
it by precessing clockwise around the impurity. Thus the 
orientation of the vortex and impurity with respect to the 
trap center greatly determines the dynamics. This com- 
bination of the interactions of the vortex with the trap 
and the vortex with the impurity then dictates that for 
a vortex to be captured by the impurity while precess- 
ing around the harmonic trap and rotating around the 
impurity, the impurity must be positioned at least closer 
to the trap center than the initial distance between the 
vortex and the trap center. 



V. CONCLUSIONS 

In summary, we studied the effects on isolated vortices 
of a localized impurity generated by a narrowly focused 
laser beam inside a parabolic potential in the context of 
Bose-Einstcin condensates (BECs). We not only exam- 



ined the dynamics (dragging and capture) of the vortex 
solutions in this setting, but also analyzed in detail the 
stationary (pinned vortex) states, their linear stability 
and the underlying bifurcation structure of the problem. 

As is already well known, the harmonic trap is respon- 
sible for the precession of the vortex around the con- 
densed cloud. We have further demonstrated that a nar- 
rowly focused blue-detuned laser beam induces a local 
attractive potential that is able to pin the vortex at var- 
ious positions within the BEC, and we investigated the 
dependence of pinning as a function of the laser beam 
parameters (width and power) for different locations in 
the condensed cloud. For a fixed beam width, we then 
explored the underlying bifurcation structure of the sta- 
tionary solutions in the parameter space of pinning po- 
sition and beam power. We found that for sufficiently 
high beam intensity it is possible to overcome the vortex 
precession and to stably pin the vortex at a desired po- 
sition inside the condensed cloud. We also studied the 
conditions for a vortex to be dragged by an adiabatically 
moving beam and concluded that for sufficiently high in- 
tensity beams and for sufficient adiabaticity it is possible 
to drag the vortex to almost any desired position within 
the BEC cloud. The possibility of vortex dragging using 
periodic, two-dimensional, optical lattices was also briefly 
investigated. Due to the lattice's saddle points between 
consecutive wells, the vortex is prone to escape to neigh- 
boring wells and, therefore, dragging with optical lattices 
is arguably less robust that its counterpart with focused 
laser beams. Finally, we presented the possibility of cap- 
turing a precessing vortex by a stationary laser beam. 
Due to the combined action of the precession about the 
harmonic trap and the precession about the localized im- 
purity, the stationary laser must be carefully positioned 
to account for both precessions so that the vortex can be 
successfully captured by the laser beam. 

This work paves the way for a considerable range of fu- 
ture studies on the topic of vortex- impurity interactions. 
Among the many interesting possibilities that can be con- 
sidered, we mention the case of more complex initial con- 
ditions, such as higher topological charge (±s) vortices, 
and that of complex dynamics induced by the effects of 
multiple laser beams. For example, in the latter set- 
ting, we might envision a situation in which a single vor- 
tex is localized to a region within a BEC by appropriate 
dynamical manipulation of multiple laser beams without 
relying on vortex pinning. Such additional studies may 
provide a more complete understanding of the physics 
of manipulating vortex arrays by optical lattices. Addi- 
tional investigations will also need to consider the role 
of finite temperature and damping, as well as the conse- 
quences of moving impurities located near the Thomas- 
Fermi radius where density is low and critical velocities 
for vortex shedding are much lower than near the BEC 
center. 

Another natural extension of our work is to study the 
manipulation of vortex lines in three-dimensional con- 
densates. It would be interesting to test whether the 
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beam could stabilize a whole vortex line (suppression of 
the so-called Kelvin modes 13411') and, moreover, change 
the orientation of a vortex [35|. Along this vein, a more 
challenging problem would be to study the pinning and 
manipulation of vortex rings by laser sheets; see e.g. [s^. 
These settings would also present the possibility of iden- 
tifying a richer and higher-dimensional bifurcation struc- 
ture. 
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